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Abstract 

We study the Single File Diffusion (SFD) of a cyclic chain of particles that cannot cross each 
other, in a thermal bath, with long ranged interactions, and arbitrary damping. We present sim- 
ulations that exhibit new behaviors specifically associated to systems of small number of particles 
and to small damping. In order to understand those results, we present an original analysis based 
on the decomposition of the particles motion in the normal modes of the chain. Our model explains 
all dynamic regimes observed in our simulations, and provides convincing estimates of the crossover 
times between those regimes. 

PACS : 05. 40. -a Fluctuation phenomena, random processes, noise, and Brownian motion 66. 10. eg 
Mass diffusion, including self-diffusion, mutual diffusion, tracer diffusion, etc. 47.57.eb Diffusion 
and aggregation 



I. INTRODUCTION 



When Brownian particles are confined along a line in a quasi one dimensional channel so 
narrow that they cannot cross each other, anomalous diffusion appears and strongly subd- 



iffusive behaviour can be observed. This p 
noticed in 1955 by Hodgkin and Keynes 



lenomenon called Single File Diffusion was first 
ll who were studying water transport through 



molecular-sized channels in biological membranes. Since then, SFD also appeared in the dif- 
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fusion of molecules in porous materials like zeolites [2H^, of charges alongpolymer chains [5[, 
of ions in electrostatic traps {g], of vortices in band superconductors [^jlsl and of colloids in 

3- 14 1 or optically generated channels [isl, Even though SFD can 



nanosized structures 



be encountered in a lot of various physical systems, most of the theoretical studies devoted 
to it are generally restricted to the simplest case : an infinite overdamped system with hard 
core interactions. 

In this paper, we present simulations results concerning finite systems of long range 
interacting particles. In particular, we focus on the dependency of the diffusion properties 
with the number of particles and the damping coefficient 7. We exhibit new behaviors, 
specifically associated to systems of small number of particles and to small damping. In 
order to interpret those results, we present an original analysis based on the decomposition 
of the particles motion in the normal modes of the chain. 

In the thermodynamic limit (infinite systems with finite density p), for overdamped dy- 



namics with hard core interactions, several analytical models |17H21| predict that at long 
times, the mean square displacement of a particle of mass m grows as FnVi with the mobility 
Fh given by 

2 [D^ _ 2 IkBT 



(1) 



with Do = kBT/{m'~f) the single particle free diffusion coefficient at temperature T, and ks 
Boltzmann's constant. If the interactions are long-ranged, only two analytical approaches 



have been undertaken so far 



22 



231], for overdamped systems in the thermodynamic limit. 



There it is proven that the m.s.d. grows as Fs\/i at long times, with a mobility Fs that 



depends on the interaction potential, through the isothermal compressibility 
spring constant K = p/kt, 



2J] or the 



where 5'(0, 0) = S{q — > 0, t = 0) is the long wavelength static structure factor of the particles. 
The diffusivity -Dg// is the effective diffusivity of a Brownian particle, taking into account 



its interactions with the other Brownian particles [25|], and differs from the single particle 
diffusivity -Do- In its last version, the expression ([2]) can be interpreted by considering that 
we can derive Fs from Fh if we replace the interparticle distance 1/p by the mean square 
displacement ksT/K of a particle in the potential well due to its neighbors. In appendix [B| 
we recover the formula ([2]) without the assumption of overdamped Langevin dynamics (E]). 

In numerical simulations and experiments, the systems are obviously finite. Periodic 
boundary conditions are used in simulations, and annular geometries in experiments. As a 
consequence of finite size effects, the asymptotic behavior at long time is always (Ax^) = 
Dj\ft. All particles in the system are then totally correlated and diffuse as a single effective 



particle of mass N x m 19|. We will recover it from our analytical model and provide 
measurements of the diffusion coefficient in good agreement with this interpretation. The 
SFD regime may nevertheless be observed in finite systems if the damping and the particles 
number are high enough, in a manner that will be precised by an analytical approach of 
finite systems dynamics. 

However, most theoretical and experimental studies have been performed for overdamped 
systems only. This is the initial assumption in the existing models for long range interactin 
particles 
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23|. The relevant experiments were generally done with solutions of colloids 
16| for which overdamping is a safe assumption. The simulations l^^] are shown in 24 1 
to be in good agreement with the theoretical prediction of Kollmann, but they also assume 
overdamping in the choice of the simulation algorithm. In order to explore underdamped 
systems, we have previously studied the diffusion in a circular channel of millimetric steel 



balls electrically charged [2J|. In this experiment, identical metallic beads are held in a 
plane horizontal condenser made of a silicon wafer and a glass plate covered with an optically 
transparent metallic layer. A constant voltage is applied to the electrodes, inducing a charge 
distribution of the beads. The condenser is fixed on loudspeakers excited with a white noise 
volta ge, and we have checked that this mechanical shaking behaves as an effective thermal 



bath 
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28 



29| . In this system the measurement of the damping constant 7 proved that the 



balls diffusion is underdamped 
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We have observed that the m.s.d. of the particles exhibit 
the SFD scaling predicted for overdamped systems, with a prefactor that is only slightly 
higher than the theoretical prediction ([2]). Unfortunately, we were not able to tune the 



damping constant experimentally. Thus, in order to investigate the specific role of damping 
on the diffusion of finite systems we have developped numerical simulations that allows easy 
changes of the damping constant. 

The paper is organised as follows : Section |IT1 is devoted to the description of the algorithm 
used in our numerical simulation. In section IIIH we present our numerical results and 
exhibit new behaviors specific to systems of small numbers of underdamped particles. We 
characterize the various regimes for the m.s.d. scaling with time, and define the crossover 
times between those regimes. In section \IV\ we give a physical interpretation of those 
regimes in the framework of our analytical model. We recover the various scaling laws for the 
m.s.d. , and provide estimates of the various crossover times, showing their dependency on 
the damping and particle numbers. We summarize our results in section |Vl Two appendices 
are devoted to complementary calculations. 

II. DESCRIPTION OF THE SIMULATION 

A. A line of particles with long-ranged interactions 

We consider point particles of mass m in a horizontal plane (xOy), submitted to a thermal 
bath at temperature T. The particles are confined by a quadratic potential in y in such a way 
that they cannot cross each other, as if they were diffusing in a narrow channel (see Fig. [T]). 
This lateral confinement is chosen to mimic as well as possible experimental situations. We 
have checked that its strength do not influence the system behavior provided the beads stay 
ordered. 



We describe the dynamics with the Langevin equation. Let = (xj, yi) be the position of 
the particle i. We do not take into account the gravity, thus describing horizontal systems. 
The particle is submitted to a confinement force —PyiGy of stiffness (5 and to the interaction 




FIG. 1. Scheme of the system 
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potential U{ri), so that Langevin equation reads 

Ti + + + —ViBy = (3) 

mm m 

with 7 the damping constant and /x a random force. In our simulations, the random force 
has the statistical properties of a white gaussian noise. Therefore, its components on both 
axes must satisfy: 

= 0, (fiyit)) = 0, = 0, (4) 

{fi,{t) fi.it')) = {fiyit)fiy{t')) = 2kBTm-i5{t - t') (5) 

where ks is Boltzmann's constant and (■) means statistical averaging. 

It is suitable to put those equations in dimensionless form, defining the following dimen- 



sionless variables: t = t/7 and x = x\/k^TJJjnrf^. It gives us 



r - + ri + VU (ri) + -^yiey = fx{t) (6) 
with the dimensionless quantities: 

and the only nonzero correlation ([5]) now reads 

= {'fiyit)lly{t')) = 25{t- V). (8) 

For the sake of simplicity, we drop the "tildes" ~ in the rest of this section. 

In order to allow a direct comparison between simulations and experiments, we take the 



same interaction potential as in our experimental set-up 2J, |28|, |29|. It reads 



f/(r.) = f/o$^A'o(^^V^), (9) 



where Kq is the modified Bessel function of second order and index 0, A and two constants. 
In principle, the sum extends to all particles, but in practice the summation is limited to 
the first five neighbours of each particle, which ensures a relative precision better than 10~^ 
and reduces the calculation time. 
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FIG. 2. Force approximation. The thick black hne represents the actual force derived from eqn. ([9]), 
the green and orange dotted lines are respectively the logarithmic and exponential approximations 
in eqn. (fTO]l . 



To decrease the computation time further, we replace the Bessel functions in the expres- 



sion of the force Ffr,- 



-VC/(ri 



Ti - r. 



ijW^i I 



by asymptotic expressions, 



X 



hx + cx \\i{x) 



-e- f 1 + ^ 
2a; V a; 




for X < 1 



for a; > 1 



(10) 



where a, 6, and c are constants. They are chosen in such a way that the force and its 
derivative are continuous for x = 1, and that the force is equal to its actual value at this 
point. Those two approximations fit very well the actual force (see Fig. [2]). 



B. Algorithm 



The simulation is based on the Gillespie algorithm |30l . |3l| that allows a consistent time 
discretization of the Langevin equation (|6]). We introduce a time step value At, which for 
consistency has to be much smaller than any other characteristic time-scale of the system. In 
dimensionless units At = 10"'^. Then the velocities Xi{t + At) and yi{t + At) are calculated 
from updating formula derived from (jS]), 



X, 



it + At) = Xiit) - [xiit) + Vt/(ri(t)) • e^] x At + V2At x /i^(t) 



Ut + A^) = Ut) 



(11) 



X 



At + V2At X fiy{t) 



where rj(t) = ^Xj(t)^ + yi{tY- The positions Xi{t + At) and yi{t + At) of all the particles 
are then calculated from 



Xi{t + M)=Xi{t)+Xi{t)M, Viit + At) = y,{t) + yi{t)At. (12) 

The components of the random noise fiy and fi^ are sampled in such a way that they have the 

properties given by equations (jlj) and ([8]), hence that they are unit normal random numbers. 
We simulate systems of particles, with periodic boundary conditions. We get from (fT2|) 
equivalent trajectories, because all beads play the same role. The system is simulated 

during a dimensionless time of 10'^, which means 10^ time-steps. The quantity of interest is 

the mean square displacement (m.s.d.) along the x direction, 

{Ax^{t)) = {[x{t + to) - x{to) - {x{t + to) - x{to))f). (13) 

where to is an arbitrary initial time. The ensemble averaging is done on every beads, since 
they all play an equivalent role. Moreover, the phenomenon is assumed to be stationary, so 
that Ax^(t) do not depend on to. For a given time t, it makes thus sense to average on the 
initial time to. Let n be the overall number of time-steps in one simulation, and Ut = t/ At. 
Then the averaging on the initial time to reads 

>>^^^^it)),), - y ^^[^"'^ + ~ (y ^ '^^^^ ~ x(^At) y ^ ^^^^ 

^ n- nt + l \ ^ n- Tit + l / ' 

where the index i is such that to = iAt, means ensemble averaging and means 
averaging on the initial time to. This way of averaging greatly improves the statistics when 
rif is smaller than n. We will use it henceforward, denoting it with the simplified notation 
(■) except in appendix lAl where it is specifically discussed. 



C. Orders of magnitude of the various parameters. 



In our simulations, we work at densities p = 33, 100 and 533 particles per meter 



temperature and interaction strengt 



and numerical simulations 
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27 



1 are such that F ranges as in experiments 
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'he 



24| 



32| . The interest of the simulations is to get acces to 



parameter values that are difficult or impossible to obtain experimentally. We vary the 
particles number between 32 and 1024. This last value is comparable to some simulations 



26l . I27j but much greater than in experiments 



24| . We vary the damping constant 



7 between 0.1 and 60 s~^, extending the experimental range toward small values of 7. 
This is to be compared to the cut-off frequency of the chain (see section IIV Ap . With our 
damping constant range, we get access to both the overdamped and underdamped dynamics 
of the particles, and are thus able to exhibit the subtle behaviors linked to underdamping. 



We simulate the same system that was experimentally studied in [2J]. In the experiments, 
the beads number vary between 12 and 37, and the density p is 477, 566 or 654 particles 
per meter. The mean interparticle distance is thus such that 1.53 < 1/p < 2.10 mm, to be 
compared to the range A = 0.48 mm of the potential. The dimensionless potential ener gy i s 



such that 6 < r < 55. The damping constant 7 ranges between 10 s~^ and 30 s~^ (see 28| . 
Fig. 6). For the experimental values of density and potential energy, the cut-off frequency 
ranges between 21 s~^ and 37 s~^. Experimentally, we are thus in the underdamped regime, 
as was already noticed in 24 1. 



III. SFD OF FINITE SYSTEMS : THE DIFFERENT REGIMES 

In this section, we present our results about the evolution of the mean square displacement 
(m.s.d.) (Ax^(t)) as a function of the time t, and focus on the effects of the particles number 

(at fixed density) and on the damping constant 7. Two typical examples are provided by 
Fig. [3l The evolution of the m.s.d. may be described by power laws (Ax2(t)) oc t", with an 
exponent a that depends on the observation time. The interpretation that will be detailed 
in part IV allow us to regroup them into three different regimes : 

• during the first regime I, < t < Tbaii, the m.s.d. grows according to i^it^; 

• during regime II, Tbaii < t < Tcou, (Ax^) may be proportional to Dt only, Fsy/i only 
or to Dt then Fs\/i, depending on the parameters of the simulation. The coefficient 
D is not necessarily the free diffusion constant Dq. When both scalings Dt and Fs^/t 
are observed, we define the crossover time r^nh between them. 

• A final regime III takes place for r^ow < (Ax^) = Dj^ft aX long times with D^ ^ D 
and Dm 7^ Dq. This final asymptotic behavior is sometimes preceded by the scaling 
H^t^ with Hn 7^ Hi, the crossover time being noted run . 
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0.1 1 10 0.1 1 10 

time (s) time (s) 



FIG. 3. Evolution of the m.s.d. (in mm^) according to the time (in s.) for a chain of 32 particles 
with density 533 m~^, temperature T = 10^^ K and interaction potential F ~ 7. The solid line 
scales as t"^, the dashed line scales as t, the dotted line scales as t^^^. 

(a) Damping constant 7 = 0.1 s~^. In this low damping case, regime II is characterized by a t 
scaling, regime III by a t'^ then a t scaling. 

(b) Damping constant 7 = 60 s~^. In this strong damping case, regime II is characterized by a \/t 
scaling, regime III by a t scaling. 



A. The small time regime (regime I) 



Regime I is defined by an evolution (Ax^) = Hit'^. It is observed in all data displayed 
in FiglU and the prefactor Hi is independent on the damping constant [see FigJ3]-a)] , on 
the system size [see Fig|l]-b)] and on the interaction potential F [see Fig|l]-c)]. In this time 
range (0 < t < Tbaii), each particle behaves independently from the others and ensures a 
ballistic flight at its thermal velocity ^Jk^T jm^ so that the constant iJi should thus be 
equal to kBT/m. The duration of this first ballistic regime is called Tbaii- From our data 
summarized in FigJU we measure the constant Hi and show in Fig. [5]-a) that it is indeed in 
perfect agreement with its predicted value. 



This behavior is obviously not observed in the simulations of the overdamped Langevin 



equation 



33| , but as already been seen in simulations of the full dynamics 3J| . In 



32( 1 , they simulate the full Langevin equation but they do not display data for a sufficiently 



small time to observe the scaling. 
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0.1 1 10 0.1 1 10 80 

time t (s) time (s) 




0.1 1 10 80 0.1 1 10 80 

Time t (s) Time t (s) 



FIG. 4. Plot of {Ax(t)'^) (in mm^) according to the time (in s.) for a density p = 533 m~^. 
Unless otherwise specified, the parameters are N = 32, T = 10^^ K, F ~ 6.8 and 7 = 10 s~^. 
Specific values are as follows : (a) 7 = 0.1, 1, 10 and 60 s^^ (blue, red, green and orange 
respectively), (b) 7 = 1 s~^ and = 4, 16, 64 and 128 (blue, orange, red and green re- 
spectively), (c) F ~ 4.4, 6.8, 9.8 and 13.4 (green, red, orange and blue respectively), (d) 
T = 10^*^, 10^^, 10^^ and 10^^ K, (green, red, orange and blue respectively). The black thick 
line is ([55]) . the dashed line is Fs^/t with the mobility Fs given by ^ and the dotted line is either 
Dj^t with D]^ given by (jlSp . in (a), (c) and (d) or Dt with D given by (|16p in (b). There are no 
free parameters in the calculations. 

B. The intermediate time regime (regime II) 

If we consider now the second regime, two different behaviors with distinct power laws 
can be observed : Fig. HJ-a) shows that for the highest values of 7, (Ax^) only grows as \/t. 
When 7 is decreased, a linear evolution in Dt appears for Tbaii < t < Tsub- For the lowest 
values 7 = 1 s~^ and 7 = 0.1 s~^, the \/t scaling completely disappears. This is clearly a 
finite size effect , as seen in Fig. HJ-b). The data displayed in this picture are recorded at a 
low value 7 = 1 s~^, and the y/t scaling is not recovered until large numbers of particles, 
typically > 128. Data from simulations with 256, 512 and 1024 particles (at constant 
density) superimpose exactly on the data for 128 particles. 

We could be tempted to explain the linear evolution in Dt by arguing that we observe 
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the diffusion of a free particle, that needs a finite time to feel the effect of confinement. If 
this should be the case, the diffusion coefficient D should be the diffusion constant for a free 
particle, which is: 



B„ = ^ (15) 

In Fig. |5]-b), we compare our numerical values of D to Dq. It is obvious that D is very 
different from Dq except at the lowest values of the density p (that is, low interactions). We 
shall see in section IIVDI that for high interactions the coefficient D actually results from 
a collective behavior of the particles. In our model [see eqn. f l42p ]. when p is high (high 
interactions), the coefficient D doesn't depend upon 7 and is given by 



D = ^,/^=^^. (16) 

In Fig.|5]-c), we see that at high density the coefficient D is actually a function of ksT/ y/mK, 
but with a numerical coefficient that is rather equal to 1/2. The dependency of D on either 
the spring constant K = U" {1 / p) or the compressibility kt = p/K indicates that collective 
phenomenons are responsible of this behavior, and that the "free particle" hypothesis does 
not account for the linear behavior observed in strongly interacting systems. The modified 
Bessel function Kq which gives the behavior of the potential U{1/ p) [see (|9])] is a very quickly 
increasing function of the density, which explains why this behavior is typical of high density 
systems. 

When the subdiffusive regime (Ax^) = F\/t is observed, as in Fig. HJ-c) and Fig. HJ-d), 
we may measure the mobility F . As seen in Fig. [5]-d), our numerical data are in excellent 
agreement with the expression Fs given in formula (|2]), even if our system is underdamped. 
We discuss in section IIVDI the relevance of this equality for finite systems. 

The numerical data displayed in Fig. HJ-c) and Fig. HJ-d) are calculated for parameters 
values that are very close (in particular, the system sizes are equal) to the relevant parameters 



of the experiments reported in 



2J|) that the value of the 



24l |. One can check (see Fig. 5 of 
m.s.d. is the same. The duration of the experiments is unsufficient to see the final D]s[t 
scaling described in the next section. 
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C. The large time regime (regime III) 



It is quite intuitive that, for very long times and finite systems, all particles become 
fully correlated and behaves as a single effective particle of mass N x m. It is a property 
of the translationally invariant mode (see section HV A|) . For small values of 7 and A^, the 
m.s.d. grows according to H^t^ , with a prefactor that is different from the constant 
Hi introduced in section IIII Al above, if at should thus be given by the resolution of the 
Langevin equation for a free particle of mass N x m : 

which is in very good agreement with our simulations as shown by Fig. [3]-e). 

At higher values of 7, and for larger systems, we only observe a linear evolution of the 
m.s.d. (Ax^) = Dj^t as shown by Fig. |4]-c) andlH-d). Similarly, should be given by the 
resolution of the Langevin equation for a free particle of mass N x m: 

This expression is in very good agreement with our results Fig[5]-f). 

The interpretation of section UV CI will show that this regime is dominated by the collective 
behavior of the particles, so that we call Tcou the time at which these collective behavior takes 
place. We recover the values of and in section HV CI from our analytical solution 
and give an estimate showing that this quadratic scaling is favored by small damping 
constants and small particles number, as is the case in Fig. |4]-a) and|4]-b). 



D. Crossover times 



Now that we know the evolution of ( Aa;(t)^) in the different regimes , we can estimate the 
different crossover times. In this section, we will proceed heuristically, defining the crossover 
times by requiring continuity of the m.s.d. for successive scalings. 

Following this method, the ballistic time Tbaii will be the time for which the curves of 
equation Hit^ and 2Dt will intersect. Depending on the expression of the diffusion coefficient 
-D, we obtain: 
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FIG. 5. a) Coefficient Hi as a function of kBT/m. 

b) Coefficient of diffusion D according to kBT/{m^) [see (fTCj) ]. The green squares represent D 
for systems of low densities (p ~ 100 and 33 m~^) and the purple circles for systems with higher 
densities (p ^ 533 m~^). 

c) Coefficient of diffusion D according to ksT ^ft^ji^^Jfrip) [the numerical coefficient is slightly 
different from that of (|16|) ]. The green squares represent D for systems of low densities (p ~ 
100 and 33 m~^) and the purple circles for systems with higher densities (p ~ 533 m~^) 

d) Mobility Fs according to 2kBT ^f^^ j ^fwmrfp [see Q]. The blue correspond to overdamped 
systems and the red triangles to underdamped systems. 

e) Coefficient according to kBT / {Nm) [see (fT7|) ]. 

f) Coefficient of diffusion according to kBT/{N^) [see ifTS]) ]. 

All axes in mm^-s~^ and logarithmic scales. All dotted lines are of slope 1. 



for weakly interacting systems. For strongly interacting systems, one has to considere the 
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effective diffusion coefficient D of (fTB]) . which gives 



ksT 2 ksT [k^ 1 iKTm 

^ball~ \ ^ball =^ T-ball~-W • (,20) 

We performed measurements of Tbaii for different values of parameters and reported them on 
Fig. E]-a) and Fig. [B]-b). One can clearly see that two different mecanisms are at stake: the 
green circles which represent systems of high densities (p ~ 533 m~^), which are associated 
to strongly interacting particles, can be easily distinguished from the red squares which 
represent systems with lower densities (p ~ 100 and 33 m~^), thus weakly interacting parti- 
cles. Formula (|T9|) seems in good agreement with the transition times of weakly interacting 
systems whereas formula fl20l) fits the values of Tbaii for strongly interacting ones. 

Let us now consider Tcou. For overdamped and large systems, it will be the time for which 
the curves of equation Fsy/t and D^t intersect, thus giving 

FsV^on = 2DNTcon =^ Tcou = = x ^e// x ^ = — • (21) 

Note that starting from equation ([2]), using the fact that -Def/ = -Do/>S'(0, 0), introducing 
the length L = N/ p of the chain, we may recast this expression to obtain 

TcoU = -75 — (22) 

^F>eff 

which is interesting as it tells us that Tcou can be seen as the time necessary for a given 
particle to diffuse over the length of the system L, with the effective diffusion coefficient 
that takes into account its interactions with the other particles. 

In the case of small damping and small systems, the SFD behavior is not observed in the 
intermediate regime [see Fig. |4]-a) and Fig. HJ-b)], being replaced by a Dt scaling with D 
given by f|T6|) . and the collective regime begins by a t"^ evolution [see IIII Cl and f|T71) ]. The 
time Tcoii may thus be estimated by 

2 ksT ksT 2 N [m 

2Dr.ou = /7^r,„, =^ ^^^^^^ = iW^^^-^i ^ ^^^^ = 7V ^^^^ 

One can see on figure |6]-c) that Tcou is in very good agreement with equations fl2T]) for 
large values of '-/N'^. For small values of '~fN'^, the data rather follow fl23p . according to the 
analysis provided in section IIV CI and particularly 
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FIG. 6. a) Measures of the transition time Tbaii according to 2/7 [Log scale; see (fT9]) ]. The green 
circles represent Tbaii for systems of high densities {p 533 m~^) and the red squares for systems 
with lower densities (p ~ 100 and 33 m~^) 



b) Measures of the transition time Tbaii according to y'Kj^m/p [Log scale; see (pOj) ]. The green 
circles represent Tbaii for systems of high densities {p ~ 533 m~^) and the red squares for systems 
with lower densities (p ~ 100 and 33 m~^) 

c) Measures of the transition time Tcoii according to TwyN'^Kx/'^P [Log scale; see (|2T]) ]. The blue 



circles represent TcoU for systems with (N/2'iT)y'^'^mKT/ p > 1 and the purple diamonds for systems 



with {N/2TT)^/-f^mKT/p < 1 (see (00])). 



d) Measures of the transition time TcoU according to N ^jKrrmjp [Log scale; the numerical co- 
efficient is slightly different from that of ()23p ]. The blue circles represent t^oW for systems 



with (A^/27r) Y^T^mK^Yp > 1 and the purple diamonds for systems with (N j^ix^^J ^'^mK'x j p < 1 
(see (liOD). 

All axes in s. All dashed lines are of slope 1. 
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IV. THEORETICAL ANALYSIS 



A. A chain of springs and point masses in a thermal bath 

In order to analyse the results presented in the last section, we have studied the Langevin 
dynamics of a chain of beads of mass m, aligned along the x axis, interacting with a pair 
potential U{x), with nearest neighbors interactions. Those two simplifying assumptions, a 
strictly ID system with nearest neighbors interactions only, are as we will see in excellent 
agreement with the actual dynamics. 

Small oscillations around the equilibrium position are described by linear springs of force 



constant K = U"{l/p) where p is the particle density at equilibrium 35|. Let x{l,t) be the 
position of particle / at time t. The equation of motion reads 

^x{l, t) = -l^Ml, t) + - [x{l + 1, t) - 2x(/, t) + x{l - 1, t)] + (24) 
aV at m m 

with the same notations as in § III A[ Let's consider a chain with periodic boundary condi- 
tions. We may introduce the discrete Fourier transform 

N N 

Xiq, t) = J2 t), x{l, ^) = ^ E t), (25) 

1=1 k=l 

with = — TT + 2Tik/N for k = 1, . . . ,N . From now on, we simplify the notations, dropping 
the dependency of the modes on the natural number k and replacing summations on k by 
summations on q. The variance of the displacement x may be calculated from the Fourier 
modes X(g,t), as (Ax2(t)) = ^^(AX2(g, t))/Ar2 ^^^j^ 

(AX2(g, t)) = ( [X(g, t) - (X(g, t))] [X(-g, t) - (X(-g, t))] > . (26) 

Let us first considere the mode g = 0, that will be noted simply X{t). It follows the 
equation 

-^X(t)+7:^X(t) = ^^^i^. (27) 

dt2 ^ ^ ' dt ^ ^ m ^ ^ 

Physically, this is the equation for a free particle of mass m in a thermal bath at temperature 

T, with damping constant 7. The solution is composed of two parts, Xd{t) which corresponds 

to the deterministic motion of the particle, and the fluctuating part X^(t) which depends 

linearly on the random forcing /i(g = Q,t). It reads 

t t' 

X{t) - X° = — [1 - e-^*l + - [df [ dt"e-^(*'-*"V(g = 0, t"), (28) 
7 ^ m J J 
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where = X(t = 0) and = X{t = 0) are the initial conditions. The corresponding 
mean square displacement contribution measured in the simulations is the double average 
defined by f|T^ . It is calculated in the appendix |X] and reads 



{AX^} 



NksT 
m'-f 



t 



-■yt 



) 



(29) 



We now considere the modes q ^ 0. Using the periodic boundary conditions x{l,t) = 
x{l + N,t), we see that each mode X{q, t) of wave number g 7^ follows the equation 

^,X{q,t)+l^Mq.t) + —(1 - cosg)X(g,t) = (30) 
aV mm m 

The roots of the characteristic polynomial associated to this equation are 



^ -2 ± - u^l^ uj] ^ 2^(1 - cosg). (31) 



Physically, the modes X{q,t) behaves as a particle of mass m in an harmonic potential 
well with pulsation Ug, forced by the random force fi{q,t). The solution of (I5U]) is readily 
obtained as 

^ Xfa.O) +._(,)jffa.O)^„^,„. ^ .,fa)Xfa.O)-X(,.0)^„_,„. ^ 

The translationally invariant mode q = scales as at small time t < I/7, then scales 
as t when t ^ I/7. The modes with nonzero wave number scales as ksTt^ /{Nm) at small 
time, and saturate toward the constant value 2kBT / {Nmuj^) at very large time. 

At intermediate time, the behavior of the modes g 7^ is determined by the relative 
values of Uq and 7/2, as illustrated in FigjTl 

- The modes such that Uq < 7/2 are overdamped. They reach their saturation value at 
a time tgat ~ The shortest saturation time is associated to g = vr, and reads 
C?" ~ l/^l if < 7/2 or t™" ~ 2/7 otherwise. 

- The modes such that Uq > 7/2 are underdamped. They oscillate at a frequency uj{q) = 

— 7^/4. Below a time which is roughly l/u{7i), all underdamped modes scale as t"^. 
Finally, all modes q ^ saturate toward 2kBT/{NmUg) at large time. 

We do not distinguish the overdamped (7 > 2uq) and the underdamped (7 < 2uq) modes, 
because the final result ( 155]) is in both cases a real function with the same formal expression. 
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n 

q 



FIG. 7. Dispersion relation oo{q) as a function of g, for q > Q (the curve is obviously symmetric 
for g < 0). The continuous line is valid for the infinite chain, the dots represent the modes for 
N = 32. The modes such that their frequency is less than 7/2 are overdamped, the other ones are 
underdamped. 

The averaging process is explained in appendix |Al One gets 



(AX2(g,t)) 



2NkBT 



1 + 



UJ+{q)e 



U}-{q)t 



(33) 



We remark that the limit g — )■ is not singular, and by taking it properly in this expression 
one recovers fl29|) . We think that it is nevertheless physically convenient to distinguish 
between the translationally invariant mode g = and the others, because they lead to 
different asymptotic behaviors. 

Using (1551) together with fl2^ . we obtain the mean square displacement as 



2kBT j t 1 . 
Nm I 7 7^ 



UJ4 



(34) 



which will be the basis for the following discussion. 

In the limit of very large damping, that is in the absence of the inertial term, a solution 
of (l24l) have been provided by L. Sjogren 



. We have thus extended his calculations to 



underdamped systems 



36| . and calculated the r.m.s. displacement in a different way to take 



into account our peculiar way of averaging (IT^ . We also extend this discussion, in the rest 
of this section, to the case of finite systems. 
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B. The ballistic regime (regime I) 



Since the small time behavior of each mode in flMjl is {ksT /Nm)t^ , and there are N 
equivalent contributions to the sum we get 

{Ax\t)) (35) 

This result is independent of A^, and thus valid in the thermodynamic limit too. Because 
of the inertial term in the Langevin equation (124|) . at very small time each particle be- 
haves independently from the others and undergoes ballistic flight at her thermal velocity 
^/k^^Jm. 

In order to discuss the duration rbaii of this regime, let us assume a finite, but large (in a 
sense to be precised later) particle number A^. The time evolution of (Ax^(t)) is determined 
by the mode g = and a summation on all modes 7^ 0. All modes in the summation flM|l . 
hence the sum itself, behaves as on a timescale such that 

t < Tbaii = min I -, ^ , ^ I . (36) 

For weakly interacting or equivalently low density systems Xbaii ~ 2/7, which was already 
heuristically derived in ( IT9|) . and observed in Fig. [6]-a). For strongly interacting or equiva- 
lently high density systems, we get Tbaii ~ I/^tt oc ^JvnjK oc \/ mKr/ p, in perfect agreement 
with our observations Fig. [B]-b) and the heuristic derivations (1TB]1 and (12U]) . 



C. The collective regime (regime III) 



This collective regime (regime III of Fig. |3l see section llll C|) is a property of finite 



systems only. This asymptotic behavior may be easily deduced from the sum 0341) . which is 
dominated by the contribution of the mode g = that scales as ^IksT / {Nni'^)] x t. This 
corresponds to the free diffusion of a particle of mass Nm in a thermal bath at temperature 
T. At very large time, the particles are completely correlated and behaves as a single particle 
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of effective mass the sum of all masses. The same result has been obtained in 

It is not difficult to estimate the time Tcow- It is the time at which the contribution of the 
mode g = dominates the sum of the contributions of all A^ — 1 other modes. Let us use 
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the simplifying Debye approximation u"^ = {K/m)q^. We thus get 

2kBT ^ 2kBT 2NkBT ^^^"^ i iys^^ 

where we have used the fact that the sum is the generahzed harmonic number -ff(Ar_i)/2,2 
which is equal to 7r^/6 up to corrections of order The fact that this time scales as 

A^^ explains why this long time regime is seldom seen in simulations (to our knowledge, the 



only exception is [3^) or experiments. This is clearly illustrated by our Fig. HJ-b), where 
we show that increasing shifts the long time regime toward longer times. This expression 
of Tcoii is equal to our previous heuristic estimate (1211) . It means that the reasoning at the 
basis of the derivation of eqn. (1571) includes in the right way the physical origin of the long 
time collective behavior of the finite chain. 

As was already quoted in section IIII C\ at very small damping constant it is possible to 
observe at large time an evolution (Ax^(t)) = H^t^. This is possible if the modes q ^ ^ 
are saturated while the mode g = still evolves as [ksT / {Nm)\t^ . This requires t < I/7 
(otherwise the mode g = scales as t), and that the contribution of the g 7^ modes to the 
sum in (15111 be less than that of the mode g = 0. Roughly speaking, we get 



\ V 7 



^^sT/iy ^ 2kBT _ NmkBT ^ V2K_ 



where we have used the Debye approximation to estimate the contribution of the q 7^ 
modes. This means that this regime is to be observed at small damping 7 and small 
particles number A^, which is precisely the case in Fig. HJ-a) and Fig. |3]-b). In this case, the 
estimate (137|) should be replaced by 

ksT , ^ 2kBT 2NkBT 1 

q^O 1 1=1 * 

This expression of Tcqu is equal to our previous heuristic estimate (125]) . allowing us to interpret 
the scaling at long time, in systems of few particles with small damping, as a collective 
behavior linked to the translationally invariant mode. This regime takes place when the 
greatest saturation time of the modes g 7^ is smaller than the time above which the mode 
g = evolves as t rather than t"^. The behavior {Ax^{t)) = Hj^t^ may thus be observable 
when 

1 1 A^ ry'^niKT , . 

<- =^ T^\P-r^<i- (40) 



7 27r y p 
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The relevance of this estimate is proved by FigEJ-c) and FigEJ-d). It also show that the 
time Tiin introduced in section [IV] is equal to I/7. 
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FIG. 8. Plot of {Ax(t)'^ > (in mm^) according to the time (in s.), in logarithmic scale, calculated 
from the analytical solution ()34p for a density p = 533 m~^, N = 32, T = 10^^ K, F ~ 6.8 and 
7 = 0.1s~^. The solid line is of slope 2, the dashed line of slope 1. This plot is to be compared 
to the relevant plot in Fig. HJ-a). We have indicated the times rtaii and Tcqu respectively given by 
and ([5U|1 . We see that both equation provides lower bounds for TcoU- 

D. The correlated regime (regime II) 

Let us now discuss the intermediate regime [s^] which, between the individual ballistic 
regime (regime I) and the collective behavior (regime III), exhibit the correlated behavior 
of the particles. While at asymptotically large time, all modes with finite (non zero) wave 
numbers have reached a constant value, in the intermediate regime, the physical behavior 
at a given time t of the chain results from a subtle balance between the modes that are 
already saturated and those that still evolve. To simplify somewhat the discussion, we will 
considere the limit 7 ^ 2u2n/N when all modes are oscillating (underdamped) and the limit 
7 ^ 2(jj^, when all modes are overdamped. 

Let us first assume a very low damping. The modes oscillates until they reach a sta- 
tionary value. The time evolution of (Ax^(t)) is due to the progressive disappearance of 
the contributions of the first oscillation of those modes q ^ that have reached their max- 
imum value. In Fig. lUJ-a) this is graphically illustrated with several underdamped modes 
(7/2 = 0.05 s~^ ^ cUjr = 25 s~^), together with the complete sum (^^. As a first approxi- 
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mat ion, 

At a given time t, the sum is dominated by the contributions of the modes that have not 
reached their first maximum, that is those modes such that t < l/uq. Let n{t) be the number 
of such modes. In the Debye approximation ujq = q^y K/m, the maximum wave number of 



those modes is Xjt^ K jm so that an estimate of n(t) is given by nit) ~ 2{N j^n^iXlt^/ K/m) 
(the factor 2 takes into accounts the modes ±|g|). The variance may thus be estimated as 

. . . (42) 



This is a normal diffusion with a diffusivity ksT/ [271 y/mK) that depends on the stiffness 
of the interaction K, showing that it is a collective effect. This is in agreement with our 
observations, see Fig. [3]-c). 

In the opposite limit of a very strong damping, all modes evolve monotonously toward 
their saturation value. As a first approximation, 

where we have used u+{q) —U-{q) ~ 7, ou+{q) ^ Ug/j and ~ —7. A mode is saturated 

at a time t > 7/w^. At a given time t, in the Debye approximation, the modes that increase 



with time are such that q < ^/m'y/ (Kt). These reasoning is graphically illustrated in Fig. |9]- 
a) where we show several overdamped modes (7/2 = 30 s^^ > = 25 s~^), together with 
the complete sum fl34p . Their number is thus n(t) ~ 2 (N/ 2Tr)^y m^y / (Kt) . The contributions 
of all such modes thus gives 



'^^bT ^ 2kBT ^ 2N jm-f 2kBT _2kBT fl^ 1/2 



(Ax'(t)) ~ =-^t X n{t) ~ =-^t X —J^ = -;^_ t'r^ = nii^ /-l^tv^, (44) 
^ ^ ^' Nm-f ^ ' Nm-f 2n \ Kt ii^/mK^ tt y mp7 ^ ' 

where we have introduced the isothermal compressibility ht in the last expression to ease 
the comparison with the exact expression of F5 ([2]). Taking into account the crudeness 
of our approximations, this estimate is extremely satisfactory because we recover the SFD 
behavior {Ax'^{t)) oc t^^^, with a prefactor that is almost the exact one. 

In the general case, the modes with wave number q such that 7 < 2uq contribute to 
a t scaling of the m.s.d. , whereas the modes such that 7 > 2u)q contribute to the SFD 
behavior, that is a t^/^ scaling of the m.s.d. The typical time r^nh at which the sub diffusive 
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FIG. 9. (Color online) Thick solid line : plot of (Ax(t)^ > (in mm^) according to the time (in s.), 
in logarithmic scale, for a density p = 533 m~^, N = 32, T = 10^^ K, T ~ 6.8, calculated from (j34|) 
for the relevant damping constant, (a) Underdamped regime, 7 = 1 s~^. The blue (gray) curves 
are, from top to bottom, the modes g = 7r/16, (7 = 7r/4, q = Tvr/IG. The dashed lines are, from left 
to right, of slopes 2 and 1. (b) Overdamped regime, 7 = 60 s~^. The blue (gray) curves are, from 
top to bottom, the modes g = 7r/16, g = '/r/4, q = 77r/16, q = Svr/S. The dashed lines are, from left 
to right, of slopes 2, 1/2 and 1. 

SFD behavior takes place is thus the inverse of this cut-off frequency, Tgub = 2/7. At a given 
particle number A^, the minimum nonzero frequency is uj2-k/n- If the damping constant 
is so small that all modes are underdamped (7 < 2uj2tt/n) 1 the m.s.d. scales as t in the 
collective regime II. Increasing the damping, at fixed A^, amounts to increase the number of 
overdamped modes, and favors the subdiffusive t^/"^ scaling for the m.s.d. This is exemplified 
by Fig.|3]-a). Increasing the particle number A^, at fixed 7, amounts to decrease the frequency 
hence to increase the number of overdamped modes, and favors the subdiffusive t^/"^ 
scaling for the m.s.d. This is exemplified by Fig. HJ-b). 

The result fl44p is only approximate, so that the numerical prefactor cannot be trusted, 
but it explains under which conditions the SFD regime may be seen in finite systems with 
periodic boundary conditions [we remind the reader that this latter is the key assumption 
leading to (l30l) ]. In appendix [HI we prove that the expression ([2]) for the mobility of long 
ranged interacting systems is valid in the underdamped case 7 < 2uJt^ too, in the thermody- 
namic limit. 
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V. CONCLUSION 



In this paper, we study the SFD of a chain of particles with long ranged interactions, 
without the simplifying assumption of overdamped dynamics. We have focused our discus- 
sion on finite size effects and the influence of low damping. We use numerical simulations 



of the Langevin equation with the Gillespie algorithm |30l . I3l| and modelize the system as 
a chain of linear springs (spring constant K) and point masses (m) in a thermal bath at 
temperature T. 

In our simulations data, we have identified several regimes for the time evolution of the 
mean square displacement (m.s.d. ) (Ax(t)^). At small times (0 < t < Tbaii), it evolves as 
(Ax(t)^) = {ksT /m)t^ . This is a ballistic flight that traces back to the inertial effects, and 
is observed whatever the damping 7 or the particle number A^. We recover this behavior in 
the thermodynamic limit (A^ — j- 00 at finite density) from our model. The prefactor of the 

scaling measured in our simulations is in excellent agreement with the theory. 

For finite systems with periodic boundary conditions, an intermediate regime (rbaii < 
t < Tcou) takes place. Depending on the respective values of the damping constant and 
the number of particles, we may observe either a diffusive behavior, a SFD behavior or 
successively both. We provide a physical explanation of those observations when we express 
the motion of the chain in terms of normal modes of oscillations. The m.s.d. (Aa;(t)^) of 
the chain results from the superposition of all those modes. The mode associated to the 
null wave number is always overdamped, and similar to the motion of a free particle in 
a thermal bath, since no restoring force is exerted on it. The modes of finite (non zero) 
wave numbers have the same dynamics as an oscillator in an harmonic well. At long time, 
the m.s.d. of all modes with non zero wave numbers saturate toward a constant value. The 
overdamped modes, which do not oscillate until they saturate, contribute to the SFD scaling 
(Aa;(t)^) oc t^^^. The underdamped modes oscillate before their saturation, and contribute 
to a linear scaling (Ax(t)^) oc t. 

In the thermodynamic limit, for systems of infinite number of particles, we exhibit ana- 
lytically the SFD behavior (Ax(t)^) = Fst^^"^ at asymptotically large time. We recover the 
mobilit y Fq that was previously calculated for long ranged interactions and overdamped dy- 



namics 
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231 ] ■ thus extending the previous calculation to systems with arbitrary damping. 



At asymptotically long time {t ^ Tcou), for a finite number of particles with periodic 
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boundary conditions, the system behaves as an effective particle of mass N x m. The 
physical origin of this behavior is the motion of the collective mode of null wave number, 
which is linked to the translational invariance of the system. For t > run > Tcou, the 
system undergoes a linear diffusion with a diffusion coefficient Dn = 2kBT/{Nm'~f). We 
show that this regime takes place whatever the value of 7. The duration of our simulations 
allows us to see this regime, and the measured diffusivity is in excellent agreement with its 
predicted value. We provide estimates of the time Tcoh ~ N'^mj/K at large damping and 
Tiin ~ 1/7 at small damping. Those estimates are in good agreement with our simulations 
data. At small particle number and small damping, a new regime takes place at times 
Tiin ^ t > Tcou. It corrcspouds to the ballistic flight of the effective particle of mass N x m, 
with {Ax{tf) = {kBT/Nm)t^. In this case, the time Tcou ~ N^mjK. 
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Appendix A: Averaging 

In this appendix, we calculate the averages used in our analysis of the simulations data. As 
explained in section Hi Bl [see eqn. ( !T4|) ]. we perform a double averaging. The first averaging 
is ensemble averaging, done on the statistical distribution of the random force and will 
be noted in this appendix (■)£ for the sake of clarity. This averaging is involved in eqns. (jl]), 
([5]) and ([8]). The second averaging is performed on the statistical distributions of and 
and will be denoted (■)o. It is obvious that those two averaging operations commute, 
and that (X^(t))o = X^{t) and {Xd{t))e = Xdif). In (l26l) . the ensemble averaging is made 
on the random force so that 

(AX2(g, t)), ^ { t) - t))^] t) - {X{-q, t)),] )^ , (Al) 



1. The mode g = 

For the translational invariant mode g = 0, all trajectories beginning at a given time 
to are equivalent so that all initial positions are equivalent. It is easy to check that 
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{{X(t) — X°)e)o = {Xd{t))Q + (X^(t))e = 0. For the same reason, the double average 
((X^(t)X^(t))e)o = because this term is hnear in and in /i(0,t). The variance is thus 



{{AX\)o = [[[Xit)-X'Y)J^ = + (X,(t)2>^ (A2) 

The variance for the deterministic part of the displacement is 

((AXj(g,t)).)o = [1 - e"-^*]^ = ^ [l - e-f , (A3) 



where the last expression is provided by the equipartition theorem for the potential energy. 

We then have to calculate the variance for the fluctuating part, {X^(t)^)^. We begin by 
calculating the following time derivative, 

^(AXj)e = ^(X,(t)X,(t)), (A4) 
= (X^(t)X^(t))e + {X^{t)X^{t)), = 2Re{X^{t)X^{t))e 

From fl28l) . we get 

t t' t 

{X,{t)X,{t)), = _l|dt'y" dt"e-^(*'-*") J dt'"e-^(*-*'")(/i(g = 0,t")f,{q = 0,t"'))e. (A5) 



The correlation for the random noise is {fi{l,t)fi{l',t'))e = g6ii'6{t — t'), with g = 2mkBT. 
Besides, /i(g = 0,t) = J2b=i f^i^^^)^ hence 

N N 

= 0, t)fi{q = 0, t'))e = ^)/^(^'' ^'))e = 2NmkBTSit - t'). (A6) 



1=1 V=l 



Performing the integrations in flASp . we get 



A(A.Y^). = 2^(l + e--'-2e-'), (A7) 



and a last integration gives the final result 



Injecting this result and ( IA3p in (IA2p gives the final result f l2^ as stated in the text. 
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2. The modes 7^ 



Physically, the dynamics of a mode X{q,t) with g 7^ is identical to the motion of 
a damped harmonic oscillator (130|) . In this case, the initial values 0) correspond to 
an initial position in a potential well and are thus not equivalent. On the other hand, 
the stationary state is quickly reached and the statistical distribution on X{q,0) is readily 
described taking different trajectories. In the data analysis, a trajectory is given by looking 
at values t + to)— X{q, Iq), and varying the time to amount to varying the initial value 
X{q, 0). The expression (126|) of the variance is thus replaced by 

{{AX^q, t)),)o = (( [X(g, t) - X{q, 0) - {{X{q, t) - X{q, 0))e)o] x 

X [X{-q, t) - X{-q, 0) - {{X{-q, t) - X{q, 0))e)o] )X ' (^9) 

As shown by (^^, X{q, t) = Xd{q, t)+X^(g, t) where ^^(g, t) is the deterministic part, linear 
in X{q,0) and X{q,0), and X^{q,t) the random part. It is not necessary to give exphcitly 
the random part, all we need to know is that it is linear in the random force fi{q,t). Since 
(X(g, 0))o = 0, {X{q, 0))o = and t)), = 0, we have ((X(g, t) - X{q, 0))e)o = 0. The 

cumbersome expression flA9p may thus be simplified to give 

{{AX'{q, t)),)o = (( [Xiq, t) - 0)] t) - X{-q, 0)] > 

= {{Xiq, t)Xi-q, + {{Xiq, 0)X(-g, O)),)^ - (AlO) 

-2Re[((X(g,t)X(-g,0)),)J 

This may be simplified further, because the stationarity implies g, t))^)^ = 

((X(g, 0)X(— g, 0))g)Q = NkBT/K. Moreover, since X(— g, 0) do not depends on fi we have 
((X(g, g, 0))g)Q = (Xrf(g, g, 0))q. The averaging is thus easily performed to give 

Appendix B: The chain of springs in the thermodynamic Umit 

In the thermodynamic limit N ^ 00, the discrete sum in flMj) may be replaced by an 
integral (taking advantage of the fact that the expressions for g 7^ are valid in the limit 
g^O), 

,^ 2. 2kBTl } ^ 1 
{Ax\t)) = / dg— 
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1 + 



u^iqjd 



u+iq)& 



-{in 



w+ (?) - (?) w+ (?) - (?) 



(Bl) 



where we used the rule (1/A^) J2q — ^ (V27r) f^^ and the invariance q — )■ 
In the hmit 7 ^ ^u^^, this integral may be expressed in closed form 



23| . This is not 



possible in the general case, but we may obtain its asymptotic behavior at large time using 



Laplace method 



. To this end, we express the time derivative 



d{Ax^it)) 2kBT 1 



dt 



dq 



m 7T J UJ+{q) — UJ^{q) 



(B2) 



For underdamped modes, u± = —7/2 ± zy^w^ — 7^/4, the large time behavior is dominated 
by exp[— 7^/2]. Since lirncug = 0, the modes with small wave numbers are always over- 
damped. The asymptotic behavior of the integral is dominated by the neighborhood of the 
maximum of uj+{q) at q = 0, with uj+{0) = 0, w+(0) = and w"(0) = —2K/{m'y). The 
leading term in (lB2p is thus 38 1 



2kBTl 



dt 



Tim 2 



.1/2 



7 



nm'jK 



211/2 



2tV2- 



(B3) 



The leading behavior of the variance, at large time, is thus 

{Ax\t))'^^-2(^^^^^^y\y'. 
To get the last expression, we have used the relation m'yD 



(B4) 

/c^T and introduced the 



isothermal compressibility of the chain, which is ht = p/K 35 1. 

We recover the behavior for single file diffusion of particles interacting with a long-ranged 
potential, as was shown by KoUmann (see 24 1 for a rewriting of Kollmann's result in 



terms of the isothermal compressibility) and Sjogren [23|. Those previous calculations were 
done under the simplifying assumption of overdamped dynamics (7 ^ 2uJt^). We nevertheless 
recover the same result because the asymptotics of integral (lB2p is dominated by the modes 
of small wavenumbers g <^ 1. Whatever the finite value of the damping constant 7, they 
are always overdamped because Uq^Q = 0. 
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